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ABSTRACT 


A simplified two-degree-of-freedom simulation in 
heave and pitch of the XR-3 Captured Air Bubble (CAB) 
testcraft is developed first ina non-linear form and 
is subsequently linearized about the steady-state 
operating point. The model is validated against the 
Six-degree-of-freedom XR-3 Loads and Motions program 
in the time domain; the linear system is then 
transformed into the frequency domain by complex 
matrix inversion. A comparison is nade between the 
linear system frequency response predicted by the 
two-degree-of-freedon model and that of the Loads and 
Motions program, as well as actual testcraft data, 
reduced to the frequency domain by a rast Fourier 
Transform (FFT) technique. This comparison of 
frequency response curves highlights a non-linear mcde 
of CAB craft behavior of possible interest to 


habitability and seaworthiness studies. 
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(ft/sec) 


Average Vertical Velocity aft of 
Center of Gravity : 
(ft/sec) 


Length of Bowseal in Contact with 
the Water Surface (Et) 
ДЕ 
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Partial derivative of vertical 
buoyant force with respect idi 
Partial derivative of vertical 
buoyant force with respect to z 
(LD AEE) 
Partial derivative Oft Vertical 


force due to planing with respect 
to Theta 
(1ЬЁ) 


Partial derivative ОЕ vertical 
ressure force with respect tO 


heta 
( lbf) 
Partial derivative of vertical 
pressure force with respect tc Mb 
( ltf/lbm) 
Partial derivative of vertical 
pressure force with respect to z 
(LEEZEE) 
Partial derivative of Bcwseal 
vertical force with respect to z 
(BEA) 
partial derivative of Bowseal 


vertical force with respect to 
Theta a 
(lbz) 
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10. 


11. 


12. 


1. 


14, 


15. 


16. 


17. 


18. 


Dpbth 


Partial derivative or Bowseal 
vertical force with respect to Mb 
(1bf/1bm) 


Partial derivative of  Sternseal 
vertical force with respect tc z 
(23 ог даа 


Partial derivative of Sternseal 
vertical force with respect ЕО 
Theta 


(128) 
Sum of Dhszi and Dhsz2 
(ПРЕРАНЕ) 
Sum of Dhsth1 and Dhsth2 
(lbf) 


Partial derivative of Plenum Air 
Mass Flow Rate with respect tc Mb 


(Sec) 


Partial derivative of Plenum Air 
Mass Flow Rate with respect to 


Theta 
( lbu/sec) 


Partial derivative of Plenum Air 
MasS Flow Rate with respect tc Z 
(lbm/ft-sec) 


Partial derivative of  Pbbar “ith 
respect to Plenum Air Mass, Mt 


2 
( lbrf/lbm-rt) 


Partial zdezivative of  Pbbar With 
respect to Theta 


2 
) / £0) 
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ШІ. Dpbz Partial derivative of  Pbbar with 
respect to 72 


3 

“ንሪ. 

20. Dpbyth Partial derivative of buoyant pitch 

ne moment with respect to Theta 
(st ihr) 
21. Dpbyz Partial derivative of buoyant pitch 
moment with respect to 2 

(405) 

Zee Dppith Partial derivative cf Planing pitch 
moment with respect to Theta = 

(t= 151) 

23. Dppmb Partial derivative ог pressure 


pitch moment with respect to Mb 
(t-il) 





24. Dppth Partial derivative of pressure 
pitch moment with respect to Theta 

(1t-1bf) 

22: DPpz Partial derivative of pressure 
pitch moment with respect to z ? 

(1р2) 

| 26. ` Dpsnb Partial derivative of Bowseal pitch 


| = 0091051۴ with respect “o ñb 2 
( гъ-Трге/ 15) 


27. Dpsthi Partial derivative of Bowseéal pitch 
moment with respect to Tneta 

(2221-1201) 

28. Dpsth2 Par а ТІ derivative of Sternseal 
pitch moment with respect to Theta 

(፲ ቲመ1 2:2) 

23. рр521 Partial derivative or Bowseal pitch 
; moment with respect to z ( 155) 
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Partial derivative of Sternseal 
pitch moment with respect to 2 198) 
р 


Partial derivative of — 
acceleration with respect to 4 


2 
( 1/(1bm-sec)) 


Partial derivative of itch 

acceleration with respect to Theta 
( sec) 

Partial derivative of pute 


acceleration witn respect to 2 р 
( 1/ (ft-sec)) 


Partial derivative ጋጀ vertical 
acceleration with respect to Mb 


2 
( ft/lbu-sec) 


Partial derivative of of vertical 
acceleration with respect to Theta 


2 
(ft/sec) 


Partial derivative Of Toren 
acceleration with respect to 2 


( sec) 
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1 ٠ INDHODUCTION 


The  f£rsquenc; domain provides a useful method of study 
As ystems such as the Captured alr Buoble Сив) # пое 
surface eiiect ship, which are excited by torces wnose 
Preece sStics of aaplituce and rre4uency day ре well-«nown, opt 
whose time domain phase relitionsnips are yuite 


Mmegeterminant. In addition, tne £reruency responses? or such 


hign-speed ships should Db» closely exanminea ог resonances 
1 


near wave encounter Írequencies ۵ lixely 


л 
= 
ct 
ጋ) 
ib 
2 
> 
ty 


experienced. u СО ICS relatively shallew crac 
ЕЕ сап сасіі/ be excited into a non-linear regien of 
Operation in which  sidewaill yaodpiny of wave/plenun upper 
emriace contact may occur; tnes2 non-linearities ay excite 
femmes ог Operation which ure of importance to naditabilicy, 
and may limit the craft's speed duriny sea-state operation 
Basten in ret [1] analyzed tre frewvuency domain resyonse 
ጠየ -1121-1 088 craft Wtilizing che Loads and Motion pro 
Шәмеіореі Dy Oceanics, Dnc., ror the bell 1003  САЗ апа 
modified extensively oy Leo and Boncal (ref I2) eni MERE 
Bene), TOES INGE дата, However, the Loads and Motion 
pregram is an extrenely ASL li cuentos 1” to. UC с 
fourier series analysis; pot меп ушар мамаз пилегіса11у 
шинээ rated alcn] the hull 2atilizing a table look-49, "once 
шэг senden rcrre@uency conponent, and rour times ior دن دوه‎ 18 
Seep, Dut whemhagher frequency components also drive tae 
variable-step incegrators to extremely small Step sizo. 
Mme fet eftect of tnis tecanlyus 15 tne з хора С 
ПЛИК асп ЕСЕП С id. a. сос  —sWMuwlaticnó ህሬ 
ЕЕ Conditions. In addition, tie time constants ot 


the  crart tesponse are relatively; lonj in бошо modes, 
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cu lariy pitch; therefore still further CPU time is 
wasted attaining steady-state conditions. This inefficiency 
of CPU time has been noted by Mitchell in ref [4], who 
documented over 40 percent of CPU time utilized py the XR-3 
Loads and Motion program to the WAVES subroutine. 
Accordingly, опе is reduced to performing a point-by-point 
frequency analysis with a single wave frequency at a time. 
Since, under these conditions, the system is approximately 
linear, and superposition is valid for sufficiently small 
wave amplitudes, the utilization of a linearized, albeit 
highly simplified, program presents itself as a most 
desirable alternative. The effect of non-linearities ray be 
added to the linear model with a describing function, and 
the system investigated in the complex plane utilizing the 
mcdified Nyquist criteria. Finally, а simplified modei 
provides a starting Point One tne Painvestigaticn of 
sensitivity of craft operation to the variation cf assorted 
parameters through the technique of root locus plotting, and 
can be the heart of the development of any automatic controi 


system, optiral or classical. 


Gerba and Thaler, in ref [5) , initiated the develcpment 
of the linearized XR-3 model in heave only. In this work, 
that model is further developed to include the pitch mode, 
first in a non-linear fora, from which a linearized version 
1s derived. The wave forcing functions acting along the 
hull of the XR-3 are developed in the form of a auasi-linear 
B vector in state-space, and the linear system 15 then 
transformed into the frequency domain by complex matrix 
inversion, with significant savings in CPU time. This 
linear frequency response is then compared with ref Г] Vo 
prcvide additional insights into the XR-3 system behavior 
not readily apparent in the more conplex Loads and Mctions 


program. 


A technique utilizing a  nmini-computer and spectral 
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analyzer is developed for the fast and efficient reduction 
of recorded XR-3 data from the time domain to tne frequency 
domain utilizing a Fast Fourier Transform (FFT). Actual 
craft data is compared with the predicted linear systen 


freguency response. 
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This chapter develops the non-linear equations of 
motions of the XR-3 in heave and pitch, following tae 
approach utilized by Gerba and Thaler in the development of 
tne simplified heave-only model in ref ЕХЕ 


А. ASSUMETIONS 





As in ref [5], the following assumptions are made: 


1. constant leakage area; 
2. equilibrium speed conditions; 
| 3. plenum pressure and volume related through an adiabatic 


process; 


The following assumptions and limiting conditions were 


removed: 


4. vertically coincident geometric center,center Of 


pressure and center of gravity; 
mee negligible pitch changes; 
6. negligible wave action; 


7. negligible added mass effects; 


and the follcwing assumptions are added: 
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8. negligible wave action within plenum; 


9. wave action due to head-on or following seas 


(negligible roll or yaw); 


In addition some simplifying assumptions regarding the craft 
gecmetry were made. Fig 1 shows the simplified gecmetry 
of the craft assumed in this model. Note that Center of 
Gravity(C.G.) and Center of Pressure (C.P.) are arbitrarily 


located and not to scale. 


В. COORDINATE SYSTEM 


Fig 1 shows the coordinate system used in tnis model. 
Specifically, the coordinate system is attacaed to the 
gecmetric center of the craft at the calm waterline, and 
oriented in a right-hand system such that x is positive 
forward, y is positive to the right, and z is positive down, 
Angles are measured such that a roll to the right, a pitch 
up and a yaw right are positive. Note, however, that 
translations in the y direction and roll and yaw are not 


included in this nodel. 
Draft, the distance of the keel pelow the calm water 


line, is the distance of the С.б. above the calm water line, 


z, plus the distance of the keel below the C.G., ZS, OTI: 


Lad =Z+2ZS 


(II-1) 


C. SIDEWALL GEOMETRY 
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The constant sidewall area as used by Gerba and Thaler 
yields a good approximation to craft motion in the heave and 
pressure modes. However ,1t is a source of strong coupling 
between draft and pitch, and a more detailed modeling is 
necessary to achieve good agreement with th2 Loads and 


Motion Prcgram, as well as known craft behavior in pitch. 


Fig 1 shows the cross-sectional detail of assumed 
sidewall geometry. The fore and aft sections are identical 
with the exception of daadrise angle, that angle being 600 
at the bow and 809 a- the stern. The actual sidewall width 
used for computation of  buoyancy forces is the effective 
sidewall width; this is the average value of the waterline 
sidewall width and the keel-line width,4s10 or Ws20 (fore or 
aft section). This is one-half the draft averaged from the 
C.G. to the bow and stern respectively, divided by the 
tangent of the appropriate deadrise angle plus the xeel-line 
width Ч510 or Ws20. This value is constrained between zero 


and Ws, the meximum sidewall width. Thus: 


Lábar 1 =Id-(L/2-Xcg) *Theta 
(112-2) 
Ws 1 =Ldbar1/(2*tan (Dr 1))+“s 10 
IF (Wsl.GT.Ws) Wsl-Ws 
ТЕ (s1.LT.0.0)Ws1=0.0 
(II-3) 


End ror the stern section: 
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Ldbar2 =Ld+ (L/2+Xcg) *Theta 


(II-4) 
WSZ =Ldbar2/(2*tan (dr2))+Ws20 
IF (Ns2.GT.Ws)Ws2-Ws 
፲፻ (፳52.5፲. 0) #852=0.,0 
(11-5) 


эн 151524 DYNAMICS 


The equations of the system dynamics are developed in 
this section and are an extension of the model of rer [5], 
modified as necessary; the remaining equations are derived 


TER first principles. 


1. Plenum Pressure and Mass Flow Rate 


Plenum pressure and plenum air mass are calculated 


uSing the adiabatic law. Using pues the draft et the bow, 
and иа аз the draft at the stern, the average draft overall 


becomes: 
Ldbar A ES. 


= (Ld- (L/2əxcg)*tan (Theta) 
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LL Ji (LZo2txXcCg)*tan(Theta)) /2 


=Ld+Xcg*tan (Theta) 


(II-6) 
which, when multiplied by the Plenum area,yields the 
submerged Plenum volume. Subtracting this fron Plenun 
nominal vclume yields the net Plenum volume, or: 

Vb =Vn-Ab* (Ld+Xcg*tan (theta) ) 
(II-7) 


The adiakatic relationship between pressure and volume 


yields: 
Ganma 
Pb =Pa* (Mb/ (Vb*Rhoa)) 
(II-8) 
and 
Pbbar =Pb-Pa 
(II-9) 


Vclume flow into the Plenum is, assuming a constant fan map: 


Gin =N* (QIO-Pbbar) 
(11-10) 
and flow cut is: 
1/2 
Qout -Cn*Al* (2*Pbbar/Rhoa) 
(II-11) 


thus the net mass flow rate within the plenum is: 
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ши (Mb) =Rhoa* (Qin-Qout) 


(11-12) 


Fig 2 shows the forces and their respective moment 
arms acting on the XR-3. The sum of vertical forces yields 


the vertical C.G. acceleration: 


2 
I (2) =(Hpres+Hbuoy+Hseal+Hdamp+Hplan+k) /M 
/ 
(11-13) 
The individual force contributions may be calculated 


separately and their total contributions summed. 


a. Pressure Forces 


The plenum pressure acting on the roof ої the 


plenum produces a pressure force acting at the C.P.: 


Hpres --Ab*Pbbar 


(11-14) 


k.e Buoyant Forces 


For purposes of calculating buoyant moments in 
the subsequent calculations, it is convenient to disscciate 


the buoyant force into two components,Hbuoy1 and Hbuoyz, the 
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bucyant forces generated by the submerged volume fore and 
aft of the C.G. Buoyant force is proportional to submerged 
volume multiplied by the density of the fluid and the 
gravitaticnal constant. This term must be doubled to account 
for both sidewalls. The average draft forward O£ the C.G. 


is: 
Табаг 1 = (La+Ld 1Y2 


=Ld- (1/2-Xcg) *tan (Theta) /2 


(11-15) 
Thus the forward buoyant force is: 
Нриоу 1 --2*Rho*G*length*draft*sidewall width 
=-2 *Rho*G* (L/2-Xcg) *Ldbari*Ws1 
(11-16) 
The aft average draft is: 
Ldbar2 =Ld+(L/2+Xcg) *tan (Theta) /2 
(II-17) 
and the aft Luoyant force is: 
Hbuoy2 --2*Rho*G* (L/2*Xcg) *Ldbar2*Ws2 
(11-43) 
and: 
Hbuoy =Hbuoy 1+Hbuoy2 
(11-19) 


C. Bow and Sternseal Forces 
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Fig 3 shows the detail of the assumed bowseal and sternseal 
gecmetry and water surface interactions. The seal is 
assumed to trail from the upper hinge point at a constant 
angle of 31.6% until it contacts the waterline. Fron that 
point on it is assumed to lie tangent to the waterline. The 
length of bowseal in contact with the water may be seen to 
be: 


Xseal1 -141/со5(31.60) 


=1.65* (Ld- (L/2-Xcg) *tan (Theta) ) 


(II-20) 


Pbbar exerts a force upon this wetted area which is resisted 
by an upward reaction from the water surface. This reaction 


is transmitted to tne craft as: 


Eseal1 =-Pbbar*Width*Xseal1 


(II-21) 


As may be seen fron Fig 3, the  sternseal 
gecmetry is basically identical to that of the bowseal; 
however, the mechanism for force generation is different in 
that the sternseal rides atop an exhausting jet or plenum 
air. Therefore, the pressure differential is not that of 
Plenum pressure over atmospheric, but over the reduced 
pressure of the high-velocity jet. This is a complicated 
function of leakage area and other factors, but may be 


assumed approximately constant. This model assumes 
2 
approximately 2.0 1bf/ft differential pressure between the 


Sternseal and plenum pressure. 


The length of sternseal in contact with tne 


water may be seen to be: 
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Xseal2 - ЖУ/Усо5(31.69) 


=1.65* (Ld* (L/2*Xcg) *tan (Theta)) 


(122) 


the sternseal force is then the differential pressure,Sealp, 


times the sternseal area ‚ог: 


Hseal2 =-Sealp*Width*Xseal2 


(11-23) 


d. Damping Forces 


Damping forces are developed for inforzation 
only, as Taylor's series expansion of these forces about the 
zero-rate equilibrium point causes these equations to vanish 


in the linear model. 
A body at motion within a fluid experiences a 
mnie tional force: 


F --C*(Rho/2)*V*|V|*S 
f f 


where C is a non-dimensional friction coefficient based upon 
f 


shape, roughness, and Reynold's Number. Consistent with the 


Six-degree-of-freedom Loads and Motion progran, — taken 
to be .9. The vertical velocity is taken to be tne vertical 


velocity of the C.G., 2401, plus the tangential velocity of 


the angular velocity, which varies along the hull as: 
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V(x) =Zdot+ (Xcg-x) *Omega 


(11-24) 


Mie” differential surface wetted area also is a function of 


X, as the depth of immersion varies along the hull, so that: 


ds =Ws* (Ld-(x-Xcg) *tan (Theta) ) 


(II-25) 


When multiplied by the absolute value of velocity, times 
itself, times the differential wetted area, the resulting 
expressicn may be integrated numerically along the null. 
However, this is most wasteful of CPU time, and the absolute 
value term renders this expression intractable to closed 
analytical sclution.. In addition, the term is guite small, 
and in fact must be, 1f pitch angles are to remain within 
the small angle assumptions, and no sidewall gapping cccur. 
Therefore, this term 1S approximated by computing the 
average velocity fore and aft of the center of gravity, and 
the average wetted surface fore and aft, previously computed 
as Wsi*Ldbar1l and  Ws2*Ldbar2. For purposes of pitch 
calculaticns to be discussed subsequently, these two forces 
are assumed to act at two-thirds the length fore and aft the 
— this yields closer agreement to the numerical 
solution by weighting the higher velocities which cccur at 


the bow and stern. Vbari and Vbar2 may now be defined: 


Ураг 1 =Zdot- (1/2-Xcg) *Omega/2 
(11-26) 
Vbar2 -Zdot*(L/2*Xcg)*Omega/2 


(11-27) 


Since there are fcur surfaces forward and aft of the center 
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of gravity the damping terms evaluate as: 


Hdamp1=-2C*Rho*Vbar1*|Vbar1]*Ldbar1* (L/2-Xcg)*Ws1 
£ 
2122-2232) 
1 е лн рат ГЭ وه د‎ #452 


(II-29) 


e. Planing Forces 


Initial time domain validation studies ccnducted 
With the model developed from the equations of this chapter, 
and the Loads and Motions program, disclosed a serious 
discrepancy between the two models in natural pitch 
frequency. Specifically, the Loads and Motions program 
predicted a natural pitch frequency of 3.2 rad/sec а 10 
kts, and 5.4 rad/sec at 20 kts, compared with a constant 
pitch frequency of 3.3 rad/sec for the linear model. This 
was taken to indicate the presence of a velccity-dependent 
source of theta feedback. Examination of ref 2] disclosed 
that Leo and Boncal had indeed utilized an added mass effect 
which had the effect of a planing force, in their 


modifications of the Loads and Motions progran. 

A body in motion in a fluid experiences a 
lifting force proportional to the square of the velocity and 
Шесту ргорсгї1опа1 to the angle of attack, ог: 

2 
Hplan азулы ы ыды 
а 


(21-00) 


where 5 is tke surface area of the body, or: 
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5 = (Ws1+Ws2) *L/2 


(21-31) 


and El is the lift-curve slope. This term is not directly 
a 


available, and is a complex function of the sidewall 


gecmetry. A value of 1.5 gave good agreement with the 
natural frequency of the Loads and Motions program, but a 
less empirical approach to this value is recommended for 


future studies. 


Ihe planing force is assumed to act at the aft 


quarter-length point.. 


The forces described above are nultiplied by their 
respective mcment arms aS shown in Fig 2 to yield the 
pitching moments about the center of gravity. Tae angular 


acceleration about the C.G. is then: 


Alfa = (2pres+Pbuoy+Pseal+Pdanmp+Pplan)/lyy 


(II-32) 


de Pressure Moment 


The pressure exerts a moment about tne center of 


Geavity, equal to: 
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Ppres --Hpres*(Xcp-Xcg) 


(II-33) 


b. Buoyant Moments 


The buoyant forces due to the immersed vclunes 
fore and aft of the center of gravity exert a pitching 
ncment alcng their moment arms, which are approximately (for 
small angles) the centroid of the rectangle from the center 


of gravity to the bow and stern: 


Pbuoy 1 ои хед) "Hbuoy1 


(II-34) 


Pbuoy2 =.5*(L/2+Xcg) *HDuoy2 


(II-35) 


c. Bow and Sternseal Moments 


The bowseal acts along a moment arm that is 
formed by the distance from the C.G. to the centroid of the 


wetted area Aseal, or: 


Pseall =-Hseal 1* (L/2-Xcg-Xseal 1/2) 


(II-36) 


The sternseal force acts on a moment arm that is 
formed by the distance from the C.G. to the centroid of tne 


sternseal wetted area, or: 
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Pseal2 =Hseal2*(L/2+Xcq+Xseal/2) 


(11557) 


4. Damping Moments 


The damping moments in pitch are the damping 
forces in heave multiplied by their respective moment arms, 
two-thirds the distance from the C.G. to the bow and stern 


respectively: 
Pdampt --(2/375(1/2-Хсауж наар 


(172358) 
Pdamp2 = (2/3) * (L/2+Xcg) *Hdamp2 


(II-39) 


€. Planing Moments 


The planing moment is the planing force acting 


at the aft quarter-length point, or: 


Eran =(L/4+Xcg) *Hplan 


(II-40) 


E. INITIAL CONDITIONS 


As in the Loads and Motion program, initial conditions 


must be carefully selected to ensure steady state conditions 
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at the beginning oí a run. Since speed does not enter into 
this program, one of the major problems of the Loads and 
Mctions program is obviated. However, the craft in this 
model is particularly sensitive to pitch, and once excited 


in this mode, is very lightly damped. 


Selection of initial conditions proceeds as follows: 


1. pitch rate and vertical velocity are defined as zero; 


Ё. a value of Pbbar is selected such that plenum air mass 


flow rate is set to zero; 


H. a value of draft is calculated such that weight and 
pressure forces are balanced by buoyant and planing 


forces; 


4. a value of Theta is calculated such that pressure 


moments are balanced by buoyant and planing moments; 


5. given the Plenum volume defined by the equilibriun 
values cf draft and Theta, and the equilibrium value of 
Pbbar, the Plenum air nass ma y be calculated 


adiakatically. 


More detail on the exact procedures for the selection of 
equilibrium conditions may be found in subseguent chapters 


dealing with the computer simulation of this model. 


43 


ዚህ 


i 





) Li u 1 


۸31۱۸ 5 
/ 


МЭТА 3015 097095715 





1 
(89801 m 211801 
ҮЗ ЇС A = -—— A E ES | 
| ВИ 7 
99 99 4j 
d ER 
93X ?|* ax ? 
------- رڅ هي‎ n — 
99Х-2/1 93፤+2፤1 


—  — —À 


— не 


зе зон ማሕል s — — — — — — — 
4 
і 
` 





~ መ — —— — — — ----? 


ASSUMED GEOMETRI Or SPHRLIRDED: XR- 3 


1 


Figure 


qu 


i 


و 





— . 0*9 es. — — — + 


11У35Н 


“መ e = s 





031 1 | 


















1" 
Yun - | 
(92Х%2/1)2 Тін | 
е 0 | 
99Х-211)-211 92Y « p 
140088... ° .. A нулан | 
| | 
Е ИРОН ЖШ د‎ 
(99X-217}-e/2 (2X«211) *£/2 
ሔ9፪- 2||1ኛ15፪- 5፤1 ИСА ЕУ СА н 
IN VOR 24114108 2 1Ү15н 


FORCES AND MOMENTS ON XR-3 


2 


Figure 


45 











X и Б 
+ (0 1с05131-8) 


BOW SEAL 





STERNSEAL 


Figure 3 - DETAIL OF BOW AND STERNSZAL 
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JIP ٢ LINERS HEAVE-=PITCH MODEL 


The non-linear equations of motion developed in the 
preceding chapter form the basis of a linearized h2ave-pitch 
model, valid for disturbances about its equilibrium point, 
defined as that point about which all rates arce zero. The 
linear model permits the use of a variety of powerful 
techniques for parametric analysis, such as root locus; it 
furthermore simplifies the design of an automatic control 
system, which could serve to further linit excursions about 


the Operating point. 


ПИШИ ТАТТОК!" 5 SERIES EXPANSION 


LOL small excursions about an operating роіп+, тһе 


non-linear relationships: 


E =F (x,y) 


q =6 (х,у) 


are approximately: 


Erd f Е (аР аў) *Dx+ (dr /dy) 7 
0 0 0 
35 dg т 


where the derivatives of F and G with respect to x and y are 


evaluated at the operating point, x and y. Cancelling the 
0 0 
equilibrium pcint from both sides and dividing through by 


time yields: 
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d/ (£) чоросу (s Ч7 (7) 
ПОЕ Жат шоо 
d =(16/4x)* д ин а T 
2 (ና) (4с/ах) Ae (x) + (16/37) vc (y) 
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i 
E 


Го 


EEJving these principles to tne  üon-linear «ode 


O 
M ceding chapter, th tg Seate 032731117  pcint i 


ሀ] 


7 


G 


e 
defined as before for the equlibrium values of th2 stat 


Variables. 


Dee PLENUS GAUGE 


rj 
23 
7 
(7 
ሀ) 
Œi 
tu 


DE RIMA دنه دل‎ 


Plenum gauge oressure has been b=en defined a 


ሀ, 


Бграг =ро-ра 


(Т1-9) 
and since Pa is a constant: 
dEkbar =áàP p 
Po, the plenum total pressure, nas been defined as: 
Gamma 
Eb (0) =21ቾ (4 (0) / (RnOa*/5 (0) ) 
Е) 


Euberetore exranáing arcund tae -quilioriun values ог чо (о) 


gung (0) : 


Ganma 
Eb (0) +dEb < Рак (Мр (0) / (апоажур (0) ) 





1 





camma 
+ (Gamma/Mb (0)) *Pa* (Mb (0) / (Rhoa*Vb (0) ) хамо 


Gamma 
- (Gamma/Vb (0)) *Pa*(Mb(0O)/(Rhoa*Vb(0)) *dVb 


(III-1) 
эг, after substituting the definition of Pb(0) and 
cancelling the equilibrium values from both sides: 

dPb = саша * (Рр (0) /“р (0) ) амр 
-Gamna* (Pb (0) /Vb (0) ) *dVb 
(III-2) 


УЕ (0) may ke further expanded, since it has been defined іп 
the preceding chapter as the empty plenum volume less tae 


net submerged volume, following the definiticn of ref [5]: 


Ур (0) =Vn-Ab* (z(0) +ZS+Xcg*Theta (0) ) 
(II-7) 
and 
dVb =-Ab*dz-Ab*Xcg*dTheta 
(TIT 3) 


Substituting the expansion of increnental volume, Equation 
(ТЇТ=3), irto the expansion for incremental pressure, 


Equation (III-2), and dividing py differential time yields: 


d Pbbar)=G * (Pb (0) /4b (0))* à Mb 
wat (Pbbar)=Ganma* (PD (0) /Mb (0) ) Sat (Мо) 
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*Gamma*Ab* (Pb (0) /Vb(0))* d/ (2) 
Jat 


*Ganma*Ab*Xcg*(pb(0)/V5(0))* G/ (Theta) 


Zar 
(сеч) 
For more compact notation tia follcwiny influence 
coefficients are defined such that: 
PUDO = Sama Ape (Pb (9) /Vo (0)):; 
(CIUS) 
оао а XJ (DD (0) ZV5(9)) ; 
(III-5) 
DUDmüUUuscauna*(Ppb(D0)/40(0)); 
11.12) 
and the pressure iif£rerantial =,uation ресолез: 
d/ (Prrar)=Dprz* 4/ (2) 
ረባ ГЕ 
рае (е ва) 
241 
зуррлож d/ (15) 
Gare 
(III-d) 


ШИГ 21111118 2155 2104 RATE ОБХІУАТІУ25 





The plenum air mass fiow rate has been defined in the 


preceding chapter as: 


d/ (Mb (0)) =Rhoa* (Qin (0) -Qout (0) ) 
sat 


(II-12) 
and 
d( d/ (МБ) )=Rhca*aQqin-Rhoa*dQout 
мат 
(III-9) 
Qin has been defined as: 
Qin (0) =N* (0i0-Pbbar (0)) 
(II-10) 
so that 
doin =-N*dPbbar 
(III-10) 
and Qout has been defined as: 
1/2 
Qout (0) =Cn*Al* (2*Pbbar (0) /Rhoa) 
(11-11) 
and 
-1/2 
dQout = (Cn*A1/Rhoa) * (2*Pbbar (0) /Rhoa) *dPbbar 
(111-11) 


Combining the two expressions for incremental volumetric 


flcw into and out of the plenum, Equations (III-10,11), with 
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the expression Есі for incremental mass flow rate, Equation 


(III-9), and dividing by differential time yields: 
2 -1/2 
d/ 2 (Mb) =~ (N¥Rhoa+t+Cn*Al* (2*Pbbar (0) /Rhoa))* 4/ (Pbkar) 
/dt fat 
(III-12) 


The influence coefficients may be defined by substituting 
Equation (III-8) into Equation(III-12) and eguating terms: 


-1/2 
Dmbz =- (N*Rhoa+Cn*Al*(2*Pbbar (0) /Rhoa)) *dpbz 
(TIES) 
- 1⁄2 
Dmbth =- (N*Rhoa+Cn*Al*(2*Pbbar (0) /Rhoa))*dpbth 
(111-14) 
-1/2 
Dmbmb =- (N*Rhoa+Cn*Al*(2*Pbbar (0) /Rnoa)) *dpbar 
(111-15) 


Thus the 11пеаг1264 differential equation for plenun mass 


flow rate may ke written: 


2 
d/ 2 (Mb) -Dmbz* d/ ( 2) 
fat /dt 
+Dnbth* d/ (Theta) 
fat 
+Dmbmb* d/ (Mb) 
/dt 


(III-16) 
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Dee HEAVE DERIVATIVES 


The heave equation has been defined as: 


2 
LOA (2) = (HprestHbuoytHsealtHdanmptHplantw) /M 
/ 


(II-13) 


Damping ісгсєѕ are neglected in the linear model.We shall 
consider the various components of heave force separately, 
and superimpose the components upon completion of all 


derivaticns. 


1. Heave Force Que To Pressure 





The Leave force due to pressure nas been defined as: 


Hpres =-Ab*Pobbar 
(11-14) 
By inspecticn, its influence coefficients become: 
Dhpz=-Ab*Dpbz 

(тр? 
Dhpt h=-Ab*Dpbth 

(ІІІ- 18) 
Dhpmb=-Ab*Dpbmb 

(III-19) 


Ed 








The kuoyant heave force equation from the preceding 


Chapter has been defined as: 


Hbuoy 


=Hbuoy 1+Hbuo y 2 


z-2*Rho*G*Ws1*(L/2-Xcg)* (Ld- (L/2-Xcg) *Theta/2) 


-2*Rho*G*Ws2*(L/2*Xcg)*(Ld* (L/2*Xcg) *Tneta/2) 


(II-15,16, 17,18) 


which may be rearranged: 


Bbuoy= 


-2*Rho*G*(Ws1l*(L/2-Xcg)*ws2* (L/2* Xcg)) * (z (0) *Zs) 


2 2 
TP*Rhno*Gf*(Wsl*(L/2—-Xcg)-Às2* (L/2t Xcg) ) *Theca/2 


and the 


Dhbz 


Dhbth 


influence coefficients may be written by inspection: 


z-2*Rho*G* (Wsl*(L/2-Xcg) *€Ws2* (L/2*Xcq)) 


(111-29) 


2 2 
=Rho*G* (Ws1* (L/2-Xcg) -Ws2* (L/2+Xcg)) 


(ЕТ1-21) 
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The expression for the bowseal heave force has been 


defined as; 


Нѕеа11 (0) --Width*Xseall(0)-*Pbbar (0) 


(11-21) 
which for small perturbations in Xseall and Pobar becomes 
dHseal 1 =-Width*Pbbar (0) *dXseal 1 


-Width*Xseali (0) *dPobar 


(TI 22) 
The expression for Xseal1 has been defined: 


Хзеа11 (0) —1.65* (z(0) tZs- (L/2-Xcg)*Theta) 


(11-20) 


which for small perturbations around z(0) and Theta(0) is 


approximately 


дХзеа1 1 -1.6 57 62-1 .6 57 )(L/2-Xcg) *dTheta 


(111-23) 


Substituting tnese relations along witn the relations for 
small perturbations in Pbbar due to z,Theta and (Мр from 
Equations (III-5,6,7} yields: 


dHseal 1 =-Nidth*(1.65*Pbbar (0) +Xseal11(0) *Dpbz) * 02 
-Width* (- 1.65* (L/2-Xcg) *Pbbar (0) 


+Xseal1 (0) *Dpbth) *dTheta 
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-Width*Xseall(0)*Dpbmb*dMb 


(T ID 24 
The influence coefficients then become: 

Dhsz1 =-Width*(1.65*Pbbar (0) +Xseali (0) *Dpbz) 

(III-25) 
Dhsth1 --Width*(-1.65*(L/2-Xcg) *Pbbar (0) 

+Xseall (0) *Dpbth) 

(111-26) 
Dhsub z-Width*Xseall(0)-*Dpbnmb 

ከ => 


The expression for the stern seal has been defined: 


Hseal2 =-Sealp*Width*Xseal2 (0) 
(11-23) 
Xseal12(0) -1.65* (Ld (0) * (L/2*Xcg) *tan (Theta (0))) 
(II-22) 
Thus for small perturbations in draft and Theta: 
dHseal2 =-1.65*sealp*Width*dz 
-1.65*(L/2+Xcg) *Sealp*Width*dTheta 
(т) 


Thus the stexrnseal influence coefficients become: 
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Dhsz2 =-1.65*Sealp*Width 


(пате 23) 


Chsth2 =-1.65* (L/2+Xcg) *Seaip* Width 


The planing force is linear in Theta , or: 


2 
Hplan SSMO ETE Sz) tL Theta 
a 


(II- 30,31) 


and the influence coefficient may be written by inspection: 


2 
Dhplth g y NUS dius 2) ED 
a 
(111-31) 


5. Sum Cf Forces 


The forces previously derived may be superiuposed so 
Bat: 


2 2 
m2 (ООРУТА а 2 (2))= 
/àt fat 


= 


(1 





(Нрге5 (0) +Нбџоу (0) +Н5еа1 (0) +Нрјап(0) ) /М 
+ (Dhpz+dhbz+Dhsz) *Kdz/M 
+ (Dhpth+Dhbth+Dhsth+Dhpltn) *dTheta/ıM 


+ (Dhpmb+Dhsmb) *dMb/m 


(III-32) 
The influence coefficients may now be defined: 
DZZ= (Dhpz+DhozZz+Dhsz) /M 
(IIT 33) 
Dzth= (Dhpth+Dhbth+Dhsth+Dhplth) / м 
(111-34) 
Dzmb= (Dhpmb+Dasmb) /M 
(III-35) 


Removing the equilibrium values and dividing through by 
differential time yield the differential equation of motion 


insthe z-direction: 


3 
а/ 3 (2) =DZZ* а/ (2) 
fat 
*Dzth* d/ (Theta) 
/dt 
*Dzub* d Mb 
Пру. еи 


(III-30) 


E. PITCHING MOMENT DERIVATIVES 
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The pitching moment equation is 


2 
a (Tbeta) -(PprestPbuoy +Pseal+Pdamp+Pplan)/Iyy 
/ 


(11-32) 


As in the previous case of heave forces, damping moments are 
neglected in the linear model; the various moments are 


considered separately and Superimposed after derivation. 


The pressure moment has been defined as 


Ppres -(Xcg-Xcp)*Hpres 


(II-33) 
ШЕ Influence coefficients may be written by inspection: 


Dppz= (Xcg-Xcp) *Dhpz 


(111-37) 
Dppih=(Xcg-Xcp) *Dhpth 

(III-38) 
Dppmb=(Xcg- Xcg) *Dhpmb 

(BE 39) 


The ruoyant moment has been defined: 
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Pbuoy =Pbuoy 1+Pbuoy2 


== (0225ЕХЄс/“Нрпоүу1/2Е(1/28Хса) “Нрцоү2/2 
2 
-Rho*g*(Wsi*(L/2-Xcg)*(Ld- (L/2-Xcq) *iheta/2) 


2 
-Ws2*(L/2*Xcg) * (Ld* (Lád* (L/2*Xcqg) *Theta/2))) 


(II- 34,35) 
and the influence coefficients may ke defined: 
2 2 
Dpbyz -Rho*G* (Wsl*(L/2-Xcqg)-Ws2*(L/2*Xcq)) 
(III-40) 
3 з, 
Dpbyth --Rho*G*(Wsi*(L/2-Xcg) *Ws2* (L/2*Xcg)) /2 
(111-41) 
3. Bow and Sternseal Moments 
The towse€al contributes a moment 
Pseal1 =-Hseal 1* (L/2-Xseal1/2-Xcg) 
(II-36) 


which for small disturtances , 1s approximately 


dEseal1 =- (1/2-Xseal1 (0) /2-Xcg) *dHseal1 
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+Нѕеа]1 (0) *аХѕеа1 1/2 


(111-42) 


EN  ተበ1100 the definition cf dXseall, Equation (III-23), 
and the influence coefficients of Hsealt, Equations 
шо -265,26,27) yields: 


дЕзеа 1 1 =(- (L/2-X*Xseal 1 (0) /2-Xcg)*Dhsz1 
+1.65*Hseal1(0) /2) *dz 
ШИЕ 2 -хсватт(0) /2-Хса) *Dhsth 1 
- 1.65 (1/2-Хса) *Hseall (0) /2) *dTheta 


- (1/2-Хзеа11 (0) /2-Хса) *Dhsmkb*ditb 


(III-43) 
The influence coefficients may now be defined: 
Dpsz1 =- (L/2-Xseal1 (0) /2-Хс4) *Dnsz1 
+1.65*Х5еа1 1 (0) /2 
(፲ ፲፲--34) 
Dpsth1 = — (L/2-Xsea1l1(0) /2-Xcg)*DhSth1]1 
- 1.6 5*(LZ2-Xcg)*Hsea11(0) /2 
(III-45) 
Dpsmb =- (1/2-Xseal1 (0) /2-Xcg) *Dnsmb 
(III-46) 


The sternseal contributes a moment: 


6 1 


~ 
`< 





Pseal2 (0) =Hseal2*(L/2+Xseal2 (0) /2+Xcg) 


(11-37) 
which for small disturbances is approximately: 
dPseal2 =(L/2+Xseal2 (0) /2+Xcg) *dHseal2 
+Hseal2 (0) *dXseal2/2 
(III-47) 
and substitution from the sternseal force influence 
S@erricients, Equations (III-29,30), yields: 
dPseal2 =((L/2+Xseal2 (0) /2+Xcg) *Dhsz2 
+1.65*Hseal2(0) /2) *dz 
+((L/2+Xseal2 (0) /2+Xcg) *Dhsth2 
*1.65*(L/2*Xcg)*Hseal2(0)/2)*dTheta 
(111-48) 


and the stern seal pitch influence coefficients may be 


defined: 
Dpsz2=(L/2+Xseal2(0) /2+Xcg) *Dhsz2 


+1.65*Hseal2 (0) /2 


(III-49) 
Dpsth2=(L/2+Xseal2 (0) /2+Xcg) *Dhsth2 


+1..05* (L/2+Xcg) *Hseai2 (0) /2 


(III-50) 
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The rlaning moment is linear in Theta: 


Рріап = (L/4+Xcg) *Hplan 
(11-40) 
and the influence coefficient is: 
Dpplth - (L/U*Xcg) *Dhplth 
(111-51) 


5. Sum Cf Moments 


The individual moment contributions may now Бе 
summed: 
Alfa(0) +dAlfa= (Ppres(0)+Pbuoy (0) +Pseal (0) +Pplan (0) ) /I yy 
+ (Dppz*tDpbyz*Dpsz) *dz/Iyy 
+ (Dppth+Dpbyth+Dpsth+Dpplth) *DTheta/lIyy 
+ (Dppnb+Dpsmb) *dMb/Iyy 
(PIE-52) 
The influence coefficients nay now be defined 
Drz= (Dppz+tDpbyz+Dpsmb) /Iyy 
(III-53) 
Dpth= (Dppth+Dpbyth+Dpsth+Dppith)/Iyy 
(III-54) 


63 





шін! 


2-27 Р7Чи6) 87266 


ri 





Dpmb= (Dppmb+Dpsmb)/I y y 


(1:11555) 


and the differential equations of motion for angular rate of 
acceleration may be written, after dividing through by 


differential time and cancelling the equilibrium values: 


3 
9/ 3 (Theta)=Dpz* d 7 
ሪ4፥ ( Jp 246 (2) 
+Dpth* d Theta 
реве yf, (meta) 
+Dpmb* d Mb 
р zm (Mb) 


(III-56) 


Е. STATE SPACE REPRESENTATION 


Fig 4 is the linear model as derived in this chapter 
presented in matrix form as: 
x =A* x 
Examination of the eigenvalues of the 9-by-9 matrix shows 
four eigenvalues at the origin (0+0j) which represent Еле 
integrated output equations for z, Theta, Mb, and Pbbar; the 
remaining eigenvalues consist of a complex pair, located in 
the vicinity of -30+30j, and a real eigenvalue in the 
ከቴ ከን 1117 ol =.5, as in ref [5], and an additional ссшр! ех 
pair near OQO+4j, representing the very lightly-damped pitch 


response. 
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The linear heave-pitch XR-3 model was validated against 
the modified Oceanics XR-3 Loads and Motions Program cf ref 
[3] prior to investigations in the frequency domain. The 


linearized program, which appears at the end of the text, 


was Simulated using the Ccntinuous Systen Modeling 
Prcgram, version III. This 1s basically a FORTRAN language 
with extremely powerful variable-step integrators and 


versatile output structure. Readers familiar with FCRTRAN 
should find little difficulty following tne CSMP TIT 
pregrams developed for this simulation; for details of this 


simulaticn language, see ref [6]. 


The Time Dcmain Program consists of three sections: 


1. Initial Conditions; this subsection identifies the 
equilibrium operating point for a given set of craft 


physical parameters. 


2. Linearization; the partial derivatives of the 
non-linear state equations derived in the preceding 
Chapter are ewaluated at this particular operating 


PORES 


3. Пупаш1сз; the weight and C.G. acceleration are 
perturbed at time zero by a step weight change, and the 
resultant linear State trajectories іп time аге 
computed. Note that the linear nodel is not 
self-starting; an initial C.G. acceleration must also 


be applied. 
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А. INITIAL CONDITIONS 


hit ial conditions involve the solution of the 
non-linear system equations of Chapter II for the 
equilibrium operating point defined at the end of the 
Chapter. 


1. Plenum Pressure 


эл» ED ‹መፎ5› mm e E Е Ёё ame — — 


A value of plenum pressure is selected so that 
plenun mass flow rate is set to zero, or, fron 
Equations (II-10,11): 


Cin =Qout 
1 
Qi0-Pbbar =(Cn*Al/N) * (2*Pbbar/Rhoa) 
(ТУ-1) 
which is a quadratic in Pbbar: 
2 | 2 2 
Prbar- (2*Q010+ (Cn*A1/N)/Rhoa) *Pbbar+Q0i0=0 
(1У- 2) 


The larger solution of this quadratic, corresponding to the 
physically unrealizable condition of reverse flow through 


the fans, is rejected. 


2. Draft and Pitch Angle 


The solution for equilibrium values or draft and 
pitch, given the previously-computed value GE Pbbar, 


requires an iterative technique, due to their close 


€7 





conning: Ап initial value of draft is first computed 
assuming zero pitch and no bow or stern seal contribution. 
The ncn-linear state equations are then solved for residual 
forces and moments due to all sources. А very smali 
perturbation in draft and pitch angle of one percent with a 
minimum default magnitude of .001 is applied and the 
differential change in residual forces and moments with 
respect to draft and Tmweta are computed. Bairstow's method 
for the zeroes of functions of two variables, as described 
in ref (7), is used to drive the residuals to zero. This 
particular froblem is somewhat slow to converge; convergence 
is enhanced by amplifying the computed changes in draft and 
pitch by 1.9, but after two hundred iterations residuals on 
the order cf ten lbf and ten ft-lbf may still be observed. 
These residual forces and moments are applied to the initial 
values of Z2dot and Aifa. 


3. Plenum Air Mass 


The equilibrium values of draft and pitca angle 


define a unique plenum volune: 


Vb -Vn-Ab* (LdeXcg*tan (Theta)) 


(IE 7) 


which, when combined with the equilibiun pressure, yields a 
unique adiakatic plenum air mass, from Equation (II-8): 
1/Ganma 
HbO -Rhoa*Vb* (Pbbar/Pa) 


(IV-3) 


The solutions to these equations define the 
equilibrium values of РЬ(0), 14(0), Theta(0), and HME (0) 
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about which the partial differential equations of the 


preceding chapter may be evaluated. 


В. LOADS AND MOTIONS PROGRAM OPERATING CONDITIONS 


The XR-3 Loads and Mction program was run at 6500 lb 
weight, C.G. 9.825 ft. ahead of the aft transom at 20 kts. 
Simulation time was 30.0 sec. for the purpose of obtaining 
reasonably guiescent values of pitch angle, draft, and 
Plenum pressure. in order to obtain identical operating 
conditions, it was necessary to adjust various constants in 
the linear model until agreement was reacned with the 
steady-state operating conditions of the Loads and Motion 


pregram. Specifically, these parameters were: 


1. seal leakage area, Al, adjusted until the two prcgrans 


agreed in Plenum pressure 
Dpyanangelitft-curve.slope coefficient, С л adjusted 
la 
until the two programs agreed in natural pitch 


frequency; 


3. minimum sidewall width, #50, adjusted until the two 


programs agreed in draft; 


4. Residual pitching moment ,Mf0, adjusted until the two 


programs agreed in steady-state pitch angle. 


The results of these trial-and-error iterations yielded 
the follcwing table of values: 
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TABLE OF OPERATING PARAMETERS 


LINFAR SYSTEM LOADS AND MOTIONS 


BSEFDOPpSDA"MSTSADY-STATEZ STEADY-STATE 


PARAMETER |VZRIABLE | VALUE VALUE 


2 2 
240.85 lbf/ft ШРОТ є 


5.02 rad/sec rad/sec 


Ws02.27 3.497 in. 


МЇ240--1875 1.21 deg 





ea COMPARISCN OF RESULTS 


The Loads and Motions program was run 3.0 secs with tne 
initial conditions as obtained above, but with a weight of 
6400 ibs. The linear model was Simultaneously run with a 
100 lb step weight reduction; data from the Loads and 
Motions program was subsequently entered into the CS#¥P III 
linear program in tabular form, for simultaneous pictting 
and error computation. For tne purposes of comparison, error 
was defined as the difference between the data from the two 
pregrams, divided by the average value of the two. This 
procedure significantly reduced the probability of 
inadvertent division by zero, as compared with normalizing 


the difference with respect to one system or the other. 
1. Plerun Pressure 


Plenum pressure and Plenum pressure error (percent) 
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plotted vs. time appear in Figs 5 and 6. The two curves are 
in quite good general agreement, with the magnitude of error 
being on the order of several parts per thousand. Тһе 
linear model exhibits a gr2ater pressure reduction in the 
initial pressure transient, than does tne Loads and Motion 
results, and this transient 1S more heavily damp2d in the 
linear model. One explanation for this differance would 
appear to be the variable leaxage area of the Loads and 
Motion program. Pitch appears to couple more strongly into 
pressure in the linear model; however, the magnitude of 
pitch is significantly greater in the linear model due to 
the absence of significant damping terms. Finally, the 
pressure in the Loads and Motion program is tending toward a 
slightly higher value, further supporting the supposition 
that leakage area has decreased in the Loads and motion 


progran. 


PAC Gap Acceleration 


A and 9 51057795020. deceleration "and с.с. 
acceleration error vs. time. They generally exhibit similar 
characteristics, but due to the larger initial pressure 
reduction of the linear model, there is a more pronounced 
overshoot in C.G. acceleration in the linear model. Tn 
addition, the stronger pitch-pressure coupling in the linear 
model also translates as stronger pitch-C.G. acceleration 
coupling. The errors appear as relatively large 
percentages; however, these are the differences of two very 


small nunmters, and should be regarded cautiously. 


3. Draft 


Draft and draft error vs. time appear in Figs 9 and 


10. As has been previously verifed by the سم‎ 
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model, draft is very well approximated by a linear mogel. 
For this particular disturbance, it appears to be slightly 
less damped and tending toward a lower steady-state value in 
the linear mcdel. Little evidence of pitch coupling is 


apparent, and errors are on the order of parts per thousand. 


4. Pitch Rate 


— — = = = — 


Fitch rate and pitch rate error are shown plotted 
vs. time in Figs 11 and 12. The linear model exhibits a 
pitch rate cf approximately twice the magnitude of the Loads 
and Moticn program; this may ostensibly be attributed to the 
Sharper pressure transient previously observed; with the 
center of pressure forward of the C.G. in both models, this 
would induce a higher initial angular acceleration in the 
linear model. In addrtıon, the linear model is 
significantly less damped in pitch; therefore this larger 
Magnitude angular acceleration would build to a larger 
angular velocity, and persist for a longer period of time in 
the linear mcdel. The frequencies compare favorably, 5.04 
rad/sec for the linear model and 5.28 for the Loads and 


Mction prcgran. 


5. Fitch Angl 


Id 


Figs 13 and 14 display pitch angle in degrees and 
pitch angle error vs. time. Both show basically identical 
characteristics; i.e., an initial pitch down induced by the 
pressure drop, and an exponential decay to a lower pitch 
angle induced by the draft change. The magnitude of pitch 
angle change reflects the larger angular velocity of the 
linear mcdel, and the difference between the damping of the 
Loads and Motion program and the linear model is quite 


dramatic. The magnitude of error, however, is on the order 
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of parts per hundred. 
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Figure 13 - PITCH ANGLE, DEGREES VS. TIME 
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V. FREQUENCY DOMAIN DEVELOPMENT 
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This chapter develops the techniques for frequency 
response of a general linear system by comolex matrix 
inversion. The ocean wave equations utilized by Booth are 
reviewed and summarized, and the forcing functions of wave 


action on the linearized xX2-3 model are developed in pitch 


and heave. Portions of Booth!s work are repeated utilizing 


the linearized model, and the results contrasted and 


compared. 


ТВК SYSTEM FREQUENCY RESPONSE 


The development of a linear system permits the 
utilization of various techniques for the development of 
frequency response characteristics which are vastly superior 
in CPU time in simulation, than the methods to which one is 
limited in non-linear models. The definition of the state 
Matrix of the previous chapter defines the state transfer 
matrix: 


-1 
H (s) < |5Т-А | *В (5) 


(7-1) 


and if the system is forced by a sinusoid, s tends toward jw 


as the transients decay and: 
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-1 


B (jw) -|jWI-A|*B(jw) 
(V- 2) 
yielding the magnitude and phase relationships: 
с (м) зд НОЮ НОМ | 
(222) 
-1 
Фи) -tan|inag(H(jw)/real(H(jw)| 
(2) 


where G(w) is the gain, in decibels, with respect to some 
arbitrary reference H(0) , and P(w) is the phase angle in 
the complex plane. Thus the frequency response of a linear 
system tc a given sinusoidal forcing function may ре 
obtained ty complex matrix inversion. While this is not in 
itself a particularly economical use of CPU time, it is 
still vastly superior to a Fourier transforn. This is 
particularly true 127 the range of frequencies 1s wide ana 
variable-step integrators driven to small step sizes which 
are not required for most of the array integration. 
Finally, it is not necessary to allow the system to attain 
Steady-state AC conditions prior to computing the frequency 
response. This is particularly important in any model cf the 


2-3, Which exhibits very light pitch damping. 


B. WAVE CHARACTERISTICS 


Booth in ref [1] provides an excellent discussion on the 
Characteristics of gravity-free surface waves with straight, 


infinitely lcng, parallel, equally-spaced wave crests, and 
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constant wave height. Such a wave advances with a velocity 
C which is cften referred to as celerity, in order to 
emphasize the fact that it is the wave front which advances, 


not the water particles. 


The regular wave may thus be modeled: 


С =G/W. 
1 
(У-5) 
Lambda -22*pi*C/w. 
i 
2 
22*pifG/w. 
i 
(7-6) 


where wis the incident wave frequency, Lambda is tne 
1. 
wavelength, and G is the gravitational. acceleration. The 


wave surface is: 


A zt) ግ... ተር ር. т) 


=A* (cos (w*x/C) *sin (w*t) +sin (w*x/C) *cos (w*t)) 
0 1. i 1 1 


(V-7) 


after substitution of Equations (V-5,6), and the approrriate 
trigonometric identities. The frequency experienced Бу the 


craft, however, is the encounter frequency: 
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я MS С– VE OS (Beta) ) ZC 
1 


2 
-wW-w*V*cos (Beta) /G 
It 


(7-8) 


where V is the craft velocity and Beta is the angle of craft 
heading with respect to the advancing wave front, Beta=0 
degrees ccrresponding to the craft and wave front travelling 


in the same direction (following seas). 
с. WAVE FORCING FUNCTIONS 


Enen aCtion introduces a vertical force in heave апд а 


ncment about the center of gravity in pitch. 
It is assumed that: 
1. wave generated forces and moments are due solely to tne 
net additional hull volume submerged by the wave; 


2. that the wave does not exist within the plenum. 


These assumptions appear justified due to the small 
magnitude of waves considered, the large inertia of the 


craft, and tte filtering action of the bowseal. 


The waveheight as a function of the x-location along the 


hull and time t may be described by: 


Wvht (x,t) -Wvht*sin(w*t*2*Pi*x/Lambda) 
e 


=Wvht* [sin («*t) *cos (x*w/C) 
e ak 
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*cos,(w*t) *sin (кем су) 
= 1. 


(215522) 


after substitution from equations (V-5) and (V-6) and tne 
appropriate trigoncmetric identities. As shown by Ра, 
the differential volume of both sidewalls immersed by the 


incident wave at time t is: 


dv (x,t) = (1/2) *Wvht (x,t) *2*W¥s*dx 


( V-10) 
and the differential neave force at time t is: 


dHbuoyw(x,t) --Rho*G*dqV (x,t) 


( V-11) 


which may, after appropriate substitution from Equations 


(У-9,10), be integrated analytically from stern to bow as: 


Hbuoyw (t) --2*Hho*G* (C/u) *vht*sin(w*L/(2*C)) *sin (w*t) 
2 1 == 


=HBUOYW*sin (w*t) 
e 


( V-12) 


The differential moment due to the differential volune 


surmerged by the incident wave at time t is: 
dPhuoyw(x,t) =Rho*G* (x-Xcg) *dV (x,t) 


=Rho*G* (x-Xcg) *Wvht (x,t) *Ws*dx 


( V-13) 


which integrates from stern to bow as: 
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1 





Pbuoyw (t) =2*Rho* WS*G*WVAt*cos (w*t) * 
e 


2 
(Со Sin (w*L/ (2*C) ) 
1 i 


* 


- (C*L/(2*w) ) *cos (w*L/ (2*C) )) 
IL 2 


*tXcg*Hbuoyw (t) 
=PBUOYW*cos (w*t) +Xcg*HBUOYW*sin (w*t) 
E 3 
( V-14) 


A particular note of caution should be introduced at this 
point with regard to the integrals for wetted area and 
wetted moments along the hull; the integrals of these 


ከ ከ ons Vanish for certain Ww . Specifically, for heave, 
1 


Equation (V-12) vanishes for all time t when: 
2 |! 

W*L/ (2*G) =n* P1 
e 


: | 1/2 
W = (2*n*P1*G/L) A LI a 


ጋ 


The pitch integral will likewise vanish for all time t if, 
foma value of n such that equation (7-12) 15 zero, ‘the 
transcendental equation obtained by substitution of equation 
(T15) into equation (V-14), and equating the latter to 
zero: 


2 1w2 2 1 
tan((n*Pi*G*L/(2*C)) )=(2*n*Pi*G/ (L*C)) 


BEL 
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is also satisfied. Thus the frequency response of the craft 
shculd exhibit a certain ripple effect as these critical 


frequencies are passed. 


The two previously-defined wave integrals HBUOYW and 
PBUOYW make up a two-dimensional B vector which couples the 
wave force onto the C.G. acceleration and pitch acceleration 


state variables, from the two forcing functions, U=sin(w*t) 
1 e 
and н“). However, it is convenient to redefine this 
= 
relationship in terms of a single forcing function since the 


Sine and the cosine are related through the derivative. In 
order to facilitate this relationship, the state matrix is 


augmented with a dummy state variable: 


= 4 U) =-w* (PBUOYW/I * Ú 
10 Wat ү ) = — 1 


( 9-17) 
Thus the new state хэн represents, after integration, the 
control ርጋ multiplied by the appropriate BE) vector terms 
of PBUOYW/Iyy. Note that if a time domain sclution of wave 


forcing action were desired, a new initial condition would 
be required, representing the initial wave buoyant nmcment; 
however, that does not enter into the solution in the 
frequency dcmain. The resulting state equations are shown 
in Fig 16. As has been previously discussed, the frequency 
response may now be computed for any frequency: 


- 1 


| | | 
G (jw) N EN | 
a е | 


) 
е 


( V-189) 
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where G(jw) is the column vector of gains for each state 
= 


variable. 


D. FREQUENCY DOMAIN SCLUTION 


The previously-defined state equations were solved in 
the frequency doma in utilizing а С5МР III computer 
Simulaticn which appears at the end of the text. This 
Simulation was based around the initialization portion of 
the time domain solution described in earlier chapters; 
however, rather than a call to integration in the Dynamic 
portion cf the model, the CSHP III reserved variable TIME 
was used cnly as a stepping variable to generate 50 


frequencies between и ва w, the specified ninimusr and 
ma X 


maximum frequencies respectively, such that: 


Delw*TIME 


(7-19) 
where Delw was defined such that when TIME =1.0, «=>. The 
e max 


incident wave frequency then becomes, for ahead seas, the 


positive sclution to the quadratic Equation (V-8), with 
cos (Beta)=-1. The influence coefficients about the 
steady-state operating point from the initiaiization portion 


were copied cnto a complex matrix Ајн, the {|{jwI-A| matrix, 
e 
and the E vector for that particular incident wave frequency 


ccrputed. The International Mathematical and Statistical 
Libraries (IMSL) subroutine ВОС. complex linear 


Simultaneous eguation solution, was used to compute the 
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vector of gains: 


-1 
G (jw) 2|jwI-A|*B(w) 
e e e 


(У - 20) 
and the resultant absolute values of G(jw) were plotted 
- 
against the encounter frequency w as decibels with respect 
= 


to the unit waveheight for heave, heave rate, and heave 
acceleration, and with respect to the waveslope, w/G, for 


ПИ, pitch rate, and pitch acceleration. Note that this 


is a linear system; accordingly, variations in waveneight do 
not affect this model's response characteristics, as was 


Ele Of the non-linear Loads and Motions program in ref |1|. 
E. COMPARISON WITH LOADS AND MOTIONS PROGRAM 


Booth in ref |1| studied the frequency response of the 
XR-3 Loads and Motion program in response to seas developed 
by the RAVES subroutine. It is of interest to compare the 
linear system developed in this study with the results of 
per 111. Duplicating his work, three runs were made 
representing ahead seas, cross seas and following seas at 
10, 20 and 30 kts and a comparison made with his results. 
Merged plots of pitch and neave were obtained ror each speed 
condition for encounter frequencies varying from 1 to 25 
rad/sec. In addition, Cal-Comp plots of pitch and heave 
were obtained as functions of the ratio of craft length to 
wavelength. ШОЕ за1 conditions апа craft parameters 


selected for initialization were identical to those selected 
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for the previous chapter with regard ue time-domain 


validation. 


TE Anead Seas Run 


ከ 11 to Fig 22 show plots of linear system 
response in heave and pitch for 10, 20 and 30 knot runs into 
ahead seas, alternating with plots reproduced from ref Гј, 
representing identical runs done by Booth with the smallest 
wave amplitudes. These were selected as representing the 
most nearly-linear mode of operation of the Loads and Motion 
progran ып sea-state. One may observe tha  аепегаї 
Similarity of the two simulations, with the speed-dependent 
resonant peak of heave in the linear model varying from an 
Encounter frequency of 3.39 rad/sec at 10 knots , to 6.45 
rad/sec at 30 kts, while the non-linear system varies from 
4.0 to approximately 5.0 rad/sec. The frequency interval 
of the non-linear study was considerably coarser than that 
utilized in the linear model, 1 rad/sec vs. less than .1 
rad/sec. The linear system is significantly less responsive 
than the Loads and Motions program, by 11 dB in heave and 6 
dB Бони ла This is attributed to the absence of 
wave-plenum chamber interactions with the plenum pressure 
Pbtar, and indicates that the presence or absence of 
significant wave action within the plenum chamber is an area 
of considerable impact on the response of the craft, and 


should be included in future studies. 


The linear system is quite sharper in its resonant 
peak in pitch, as should be expected from the significantly 
reduced pitch damping previously demonstrated by the linear 
model in the time domain validation. This accounts in part 
for the tendency of the linear model to decrease in response 
slightly from an encounter frequency of 1 rad/sec, until 


just before the onset of resonance, unlike the Loads and 
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Motions data; the half-power points in the Loads and Motions 
frequency resonance are more widely separated than are those 
ШЕ the linear system, and this in part also accounts for 


this difference. 


The linear system shows the distinctly periodic 


frequency response predicted, with a sharp null occurring at 


incident wave frequencies, W ӘР Mg Uu апа 5938 


( 
1 
rad/sec. These frequencies correspond to a value of n, in 


TST 15), сі 1, 2 апа 3. These nulls are quite sharp, and 
are not as well defined in the Loads and Motions studies due 
to the coarser frequency resolution, but are nonetheless 


mest ple in all runs. 


The Loads and Motions studies also demonstrate an 
interesting form of non-linear response with regard to 
resonant pitch frequency. For small wave amplitudes, less 
than .1 ІС, the Loads and Motions resonant pitch frequency 
varies with craft speed in nearly an identical manner with 
the linear system, with the slight exception of the 30 kt 
study; however, the only available data for that run was at 
a wave amplitude of .2 ft. For wave amplitudes greater than 
„1 ft, the resonant peak occurs at a frequency of U rad/sec, 
regardless of craft speed. A tentative explanation for this 
behavior 15 herein advanced: as the craft resonant 
frequency is approached, bowseal gapping commences for the 
larger wave amplitudes, with attendant loss of plenum 
airmass; the loss of seal force and plenum pressure cause 
mre craft to settle back into the water, but a fixed amount 
of time must elapse before the fans can replenish the plenum 
airmass and restore normal system operation. Thus the 4 
rad/sec period observed ror larger wave amplitudes 
represents the fastest tine (.66 sec) over which gapping 


can occur and the fans sufficiently replenish the airmass to 
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permit it to recur. This non-linear response 15 somewhat 
analogous tc a relaxation oscillator, and represents an 
interesting area fcr future study, as such a non-linearity 
night easily be modeled in the Nyquist plane, utilizing the 
modified  Nyguist stability criteria to predict this 


non-linear behavior. 
a. 10 Kts Ahead Seas 


Figs 17 and 18 show the superimposed heave and 
pitch response of the linear and Loads and Motions progran, 
respectively, for an ahead-seas run at a craft velocity of 
10 kts. The linear model represents the response of the 
craft in heave and pitch at 50 discrete frequencies ranging 
exponentially from 1.0 to 25.0 rad/sec. The non-linear data 
reproduced from ref (1) represents 8 discrete frequencies 
over that same range at 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 6.05 
and 10.0 radysec, at a waveheight of .2 it. Although this 
waveheight was the smallest used by Booth for runs at these 
speeds, it is still sufficiently large to excite the 
aforementioned non-linear response. This behavior makes 
itself manifest as the 4.0 rad/sec resonance in pitch, 
ccmpared with 3.39 rad/sec for the linear model. The Loads 
and Moticns program shows, for small weight disturbances at 
these speeds, a natural pitch period of 1.64 secs, or 3.7 
rad/sec. As previously noted, the magnitude of the linear 
System pitch response agrees fairly weli, within 3 dE, о: 
the non-linear Loads and Motions data, with a sharper peak 
due to the lighter dinear damping; the heave magnitudes 
differ considerably, by -15 dB, due to the absence of flenun 


waves in the linear systen. 
with the exception of the heave magnitude 


difference and the shift in resonance due to the  ncn-linear 


behavior, tte shape of the heave response curves of the two 
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models ccmpare quite favorably. А Shallow na of 
undetermined origin, appears in both heave curves at 3.0 
radysec, just slightly before the onset of resonance. A 
very sharp null occurs at 8.36 rad/sec in bcth models; this 
is at an incident wave frequency of 3.15 rad/sec, or a 
wavelength of 20 ft. AS predicted by Equations (V-15,16), 
further nulls occur at a wave incident frequency of 4.46 and 
5.58 rad/sec, corresponding to wavelengths of 10 and 6.6 
feet respectively, but the Loads and Motion  enccunter 


frequencies did not range sufficiently high to observe then. 


Pitch also shows general agreement; but  Bcoth's 
Loads and Motion data shows a null at 6.0 rad/sec, vs. 11.54 
rad/sec for the linear system. This null may, however, be 
related to the non-linear response previously described, 
representing a freguency at which the phase relaticnship 
between plenum airmass replenishment and pitch buoyancy due 
to waves is such that cancellation occurs. те 15 
suspiciously located at а 3/2 multiple of the non-linear 
response frequency of 4.0 rad/sec. In addition, Booth noted 
this identical discrepancy in comparison with the linearized 
Aercjet XR-3 model; the Aerojet model also predicted a null 
in the vicinity of 12 rad/sec, which Booth was unable to 
observe with the Loads and Motions program, but which is 


distinctly visible in the linear model at 11.54 rad/sec. 


bz Z0 Kts Ahead Seas 


Figs 19 and 20 show the linear and Loads and 
Motions heave and pitch frequency response for a 20 xt run 
into ahead seas. The Loads and Motions data is for a wave 
height of .1 ft and is more nearly linear in behavior than 
previously cbserved at 10 kts. The range of frequencies 
Spanned by the linear system is the same as before; the 


Loads and Motion data includes two additional points at 15 
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and 20 rad/sec in heave response. The linear pitch response 
is smaller than that observed at 10 kts, being abcut 6 dB 
less than the non-linear data, while the heave response 
remainS consistently less than the non-linear data by 12 dB 
as at 10 kts. 


The shape of the two sets of frequency response 
data is in guite good agreement for these conditions, as the 
Loads and Motions is apparently in a nearly-linear ucde of 
operation. The pitch data shows a distinct flattening at 
the peak of resonance in both models, although this cccurs 
somewhat «earlier in the non-linear data (4.0 rad/sec) than 
in the linear дата (4.68 rad/sec), due to the more 
sharply-peaked pitch response of the latter. Both 5 
resonate at 5.0 rad/sec, and fall off above resonance with 
identical rcll-offs of 12 dB per octave, after the initially 
steeper roll-off of the linear system from peak resonance. 
Heave shows the characteristic null in both models at ап 
encounter frequency of 14 rad/sec; this is a wave frequency 
of 3.15 rad/sec, or a wavelength of 20 ft. The linear 
system predicts a pitch null at 20.0 rad/sec, but tne Loads 
and Motions data does not extend sufficiently high to 


observe this. 


C. 30 Kts Ahead Seas 


Figs 21 and 22 show the linear and ncn-linear 
frequency response for a 30 kt ahead-seas run. As in the 10 
kt run, the smallest wave amplitude utilized in Bcoth's 
study was .2 ft; however, it is considered that tne 
nose-down planing moment which is quite strong at high 
Speeds restrains bowseal gapping, and the Loads and 5 
data represents conditions for the non-linear response which 
is just cn the verge of becoming established. This is 


supported by the observation that pitch resonance cccurs 
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between the frequency previously observed as the 
characteristic non-linear response, and linear resonance, 
between 4.0 and 5.0 rad/sec; furthermore, the heave 
resonance is weakly peaked at a slightly higher frequency 
than pitch response, approximately 5.0 rad/sec, still closer 
to the predicted linear resonance at 6.46 rad/sec. 


The linear modei shows still lower magnitudes of 
pitch response, 6 to 8 dB below the Loads and Motion data. 
The heave magnitudes remain consistently 12 dB below the 
Loads and Motion data. 


Fitch resonates at 6.46 rad/sec in the linear 
system, and between 4.0 and 5.0 rad/sec in the non-linear 
Loads and Motions data.  Roll-off in pitch remains the sane 
above resonance, approximately 12 dB per octave in both 
mcdels, but heave response in both models falls cff mors 
shallowly, approximately 10 dB per octave in both ncdels. 
The heave null due to the 20 ft wavelength is still apparent 
in the linear model, representing an encounter frequency of 
20.0 rad/sec, but the Loads and Motions data did not extend 


sufficiently high in frequency for this null to be observed. 
d. Frequency Response vs. Wavelength 


Figs 23 and 24 are plots of heave response, vs. 
the ratio cf craft length to wavelength for the linear 
system, and vs. the ratio of wavelength to craft length for 
the Loads and Motion data reproduced from ref [1]. Both 
plots are superimposed with the data from all three velocity 
mens of ТО, 20, and 30 kts. The difference in scale 
selected for plotting in the linearized model precludes 
crcwding at the low end of the x-axis, since the bulk cf the 
frequencies studied, 3.15 rad/sec and above, are shorter 
than craft length. 
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The linear system heave response shows distinct 
nulls for waves which are integer multiples or submultiples 
of the craft length, and пау be observed at L/Lambda= 
174,1.0,2.0, and 3.0; less distinct nodes may be observed at 
L/Lambda= 3/2, 5/2, апа 7/2. The corresponding nulls at 
Lambda/L=4.0,1.0, and .5 are visible in Fig 24, but scale 
crowding and the coarser frequency resolution utilized by 
Booth render these less obvious. The sharp node at 
L/Lambda=.3 corresponding to a wave frequency of 1.75 
rad/sec is nct the result of the periodic hull inteorals; 
this wave frequency corresponds to the enccunter frequency 
ot 3.4 to 6.6 rad/sec, and is the craft's natural pitch 
frequency. It may te observed, also as a rode, in the pitch 
response. The natural pitch frequency variation with speed 
is obscured in this plot, since as the speed increases, the 
natural pitch frequency increases; but this is an encounter 
frequency, and the wave incident frequency which prcduces 
this enccunter frequency decreases with increasing speed. 


Therefore, the two effects nearly cancel. 


Figs 25 and 26 show pitch response for the 
linear system and the non-linear loads and Motion data of 
ref (ij, Fictted as was heave. Note that the nulls and 
nodes are roughly 180° out of phase, a node in pitch 


corresponding to null in heave and vice-versa. 


2. Сто55-5еа5 Runs 


a — ጨ= 4127 v cm‏ سه سی کس ہے سد 


Figs 27 to 31 are presentations of heave and pitch 
frequency response of the linear system and the Loads and 


Motions data from ref ከ] at 20 and 30 kts under  cross-seas 


conditions (нан). Comments made regarding the the 
e 1 


ahead-seas runs are generally applicable here, but tne 
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comparison of the two systems is in no manner nearly as 
striking as that for the ahead-seas runs. This is 
attributed to the significant effects of roll and yaw nodes 
on the craft dynamics under these conditions, which are 


completely neglected in this linear model. 
3. Following Sea 


Booth acknowledged little agreement between his work 
with the Loads and Motion program and the Aerojet XR-3 


linear model frequency response curves; this is apparently 


due to an error in calculating the encounter frequency. 
Higher frequency waves propagate at a slower velocity, so 


that for a craft traveling at a given velocity, as the 


incident wave frequency w increases, a frequency will be 
1 


attained above which the craft will actually outrun the 


following wave, and plow into such a wave bows-on. Thus the 


craft will actually experience an ahead-sea condition above 


that wave frequency. The w corresponding to the maximum 
i 


positive w that the craft may experience for following seas 
e 


~ 


(cos(Beta)=1.0) is given by: 
2 
d/dw(w-V*w/G) =0.0 
ТӨРІ і 


(7-21) 
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м =G/ (2*V) 


(7-22) 
Кес си запад Equation (V-22) into (V-8) yields: 
2 
м ከጋሪ ሚን ИС ДУДА БРАНИ) 
emax 
=G/ (4*V) 
(V-23) 
incident wave frequencies above or below the w corresponding 
1 
to w result in a lower w, and if wis increased 
emax - È 


sufficiently, past w=G/V, w will become negative, indicating 
1 e 


m the crart is overrunning the wave. Accordingly, a great 
deal of judiciousness must be utilized when interpreting the 
results of a following seas run, when plotted against 


encounter frequency w. 
e 


Accordingly, the following-seas frequency spectrum 
was generated in a manner different from that utilized on 


previous runs. W was varied exponentialiy over the range of 
1. 


1 to 25 rad/sec, and the corresponding w was calculated. 


` 


Hence, the results are not comparable directly with the 


posults of Booth, but are presented in Fig 32 to Fig 33 for 


information only, at v=20 kts. 
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Figure 17 - LINEAR SYSTEM FREQUENCY RESPONSE, 


10 KTS, AHEAD SEAS 
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Figure 20 - LOADS AND MOTION FREQUENCY RESPONSE, 
20 KTS, AHEAD SEAS 
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This chapter develops the techniques for fast efficient 
reduction of XR-3 testcraft data from the time domain to the 
frequency dcmain. Some more prominent noise sources аге 
identified, and the mode of noise transmission is 
hypothesized. A very limited amount of pitch data is 
collected and closely examined using these techniques and 
compared with the theoretical development of the preceding 


Chapters. 


A. DATA ACQUISITION SYSTEM 


The data acquisiricn system consists of the standard 
semi-permanently installed  Pemco tape recorder. This is a 
14-channel FM tape system, with output 0-1 volt positive. A 
19-channel ¡patch panel in the tape recorder compartment, 
located just aft of the cockpit, permits selection of the 
data to be recorded;voice edge track recording directly from 
the craft's radio/ICS system permits relatively simple 


documentation of experimental runs to be maintained. 


В. FREQUENCY DOMAIN DATA REDUCTICN SYSTEM 


The Eio-Engineering Laboratory, Electrical Engin¢ering 
Department, NPS, Monterey, CA, has a very powerful systen 
for spectral analysis of analog data. This systen, 


consisting СЕ а PDP-11 Mini-computer and  Time-Data 
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Corporation Superpanel-50 Spectrum Analyzer, was utilized 
for time-to-frequency domain data reduction of the analog 
tape data. This system is capable of performing a Fast 
Fourier Transform (FFT) on analog data of magnitudes of +5V, 
frcm a maximum upper frequency of 10 Hz to 10 MHz, for up to 
8094 discrete time samples. However, caution must be 
exercised tc avoid memory overflow and/or program interrupt 
of the EDP-11 for samples in excess of 2048. In addition, 
this entire process may be repeated up to 2048 times, and 
the resultant Fourier coefficients averaged, for further 
reduction cf transient frequency components. The resultant 
transformed data may be output as Fourier coefficients, 
power spectral coefficients, and as an amplitude histcgran; 
output devices include a storage CRT and a Hewlett-Packard 
7004-B X-Y plotter which will reproduce the display on the 
CRT con ccomand from the panel. In addition the time domain 
data being analyzed may be observed on a time-base storage 
oscillosccpe, and monitored throughout the transformation 
precess. Processing times vary, of course, with upper 
frequency, rumber of samples, and number of averaging 
"windows", from less than a second to several minutes. FOr 
further  infcrmation see ref [в], the operating instructions 
for the Superpanel-50. Fig 34 shows the details of the 


setup for data reduction to the frequency domain. 


С. FREQUENCY DOMAIN ANALYSIS 


XR-3 cata was codlected on an experimental run conducted 
20 January, 1977 at Lake San Antonio for the purpcse of 
performing frequency domain analysis, using the 
previously-described FFT techniques. Two analyses were 


performed: 


1. a broad-band noise analysis, examining the freguency 
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Ера Ка channels up to 200 Hz, as various 
suspected noise sources were brought up and shut down, 
or varied in power setting; the purpose of this 
experiment was to identify affending noise sources by 


their characteristic signature in the frequency domain; 


2. a lcw-frequency analysis of craft characteristics in 
pitch and pitch rate, for 15 kt runs into, across, and 
with the prevailing wave directions. The frequency 
spectra of these channels were analyzed from .025 to 
25 82, for the purpose cf comparison with reí [J] and 


the werk cf the preceding chapters. 
nd Noiss Analysis 


Ncise data from the XR-3 was analyzed over a 
frequency range up to 200 Hz and 512 samples, averaged over 
8 windows; this gave a total sample time of approximately 30 
seconds. A Butterworth filter was used to reject signals 
abcve 200 Hz to prevent aliasing of higher frequency 
components into the lower range of interest due to 
discretization. In addition, a battery/rkeostat combination 
was used tc bias out the 1/2 volt offset of the FM tape, 
Since , as previously noted, the FM tape output range is 0-1 
vclt positive, while the input to the Analog-to-Digital 
converter of the Superranel-50 is +5V; this considerably 


improves amplitude resolution of the transformed data. 


A representative output of frequency spectra appears 
in Fig 35; the remainder may be examined in Appendix A. The 
data from all channels are basically identical in nature, 
and too repetitive to include in the main text. AS may be 
expected, the craft main engines and the 60-Hz power supply 
contribute tte major share of noise. AS power is varied, the 


main engine noise may be observed to vary from 20 Hz to 100 


123 


u 





Hz, while the power supply contributes significant noise at 
55-60 Hz and its harmonics. Following these sources in 
order of magnitude are the lift fans, which may be observed 


Enocthe vicinity of 40 Hz. 


An attempt was made to ascertain the mode of 
transmission of this noise, specifically, whether by 
electrical means, or by mechanical vibration of the sensor 
package. Accordingly, the gyro package containing pitch, 
roll, and yaw, and pitch rate, roll rate, and yaw rate gyros 
was removed and placed on a Ling shaker table, and exposed 
to vibrations of 1/2 to 1G, over a frequency range of 5-200 
Hz. The RMS shaker table displacement was measured with a 
calibrated Bentley accelerometer, and compared with the RMS 
voltage output of the individual gyros. It was hoped to 
Obtain frequency response curves of the gyro packages, but 
they seemed quite insensitive to vibrations in this 
frequency range, and the frequency curves showed significant 
scatter due to the low signal-to-noise ratio. They were not 


considered worthy of inclusion. 


It was concluded, on tne basis of this experiment, 
that mechanical vibration соза not transmit such 
overwhelming amounts of noise aS was observed being injected 
into the XR-3 data collection system. The mediun of 
transmission is thus electrical; this concluSion is further 
supported by the observation that the eight engines 
installed aboard the Хана have no ignition noise 
Suppressors; likewise, the 60 Hz power supply is unfiltered. 
Furthermore, while bonding straps run the length of the 
craft, the possibility of extensive AC ground loovs is 
always present in a non-conducting hull, such as the XR-3. 
These areas should rec2ive some close scrutiny in the 
future, for a tenfold reduction in noise injection into the 
data collection system should be possible, by means of 


ignition noise suppression, and the installation of .1 uF 
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Sdpdcimercmedt each AC outlet, at a relatively insignificant 


cost in time and money. 


2. Low-Frequency Craft Dynegic 


» 


Three runs were made during the aforementioned data 
collection runs, in which pitch and pitch rate were recorded 
for 15 kt runs into, across, and with the prevailing seas. 
Qualitative sea conditions at the time of observation were 
noted as wave heights of 3 to 6 inches, with a wavelenoth of 
approximately 6 feet; this wavelength represents an incident 
wave frequency of approximately 2.2 rad/sec, corresponding 
to encounter frequencies of 6.2, 2.2, and -1.8 rad/sec, гог 
ahead, cross, and following seas at 315 kts. Data was 
reduced in a manner similar to that used for wideband noise, 
with the exception that a maximum frequency of 25 Hz, 1024 
samples with no averaging was utilized. i А difficulty 
developed in that a Butterworth filter with cut-off below 
100 HZ was not readily available at the time. However, when 
спе is constrained from lowering the bridge, it ma y 
nevertheless be possible to raise the river. In this 
particular case, the tape was run at four times recorded 
speed, or 7 1/2 inches per second. The FFT was set for 100 
Hz, as was the Butterworth filter, ana tne desiredc 
resolution of 1024 samples selected. This procedure nad the 
beneficial side-effect of reducing the real elapsed 
prccessing time of the Superpanel-50/PDP-11 combination by a 
factor of four. The resultant transformed data was 
scale-expanded after computation, so that only the portion 
веси 0—5 Н2 "(0-1.25» асъшай Hz) was plotted. 


These admittedly unorthodox and awkward procedures 
were necessitated, not only by the lack of Suitable 
filtering, but also by difficulties experienced in obtaining 


prcper operation of the PDP- 11 for very low frequencies at 
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that particular time. Long sample times combined with 
frequent inexplicable program interrupts rendered this 
particular sessicn most frustrating. However, ventilation 
and overheating problems were apparently the major source of 
difficulties experienced with this particular exercise, and 
future data reductions in the frequency domain should 


require ncne of these awkward circumlocutions. 


a. Pitch and Pitch Rate-Ahead Seas Run 


Figs 36 and 37, the reduced FFT of pitch and 
pitch rate for the 15 kt ahead-seas run, шау be compared 
with Figs 38 and 39, the linear system theoretical response 
in pitch and pitch rate. The linear system was plotted as 
dB with respect to waveslope logarithmically against 
encounter frequency, on a range of .1 to 7.85 rad/sec. The 
observed data are relative magnitudes plotted linearly vs. 
encounter frequency, on a range of 0 to 7.85 rad/sec. The 
unusual scale was necessitated for the observed data, to 
obtain the conversion to rad/sec; in Hz, the scale is .125 
НС Ln. A qualitative comparison of the two plots shows a 
rather remarkable similarity; the FFT of pitch shcus a 
falling off in response from the very low frequency end to a 
local minima at 2.0 rad/sec with the onset of  rescnance 
occurring in the vicinity of 2.75 rad/sec; resonance peaks 
at 3.9 radysec and is completed at approximately 5.0 
rad/sec. Two small nodes may be observed at 2.36 and 6.23 
rad/sec, the latter corresponding to the enccunter frequency 
of the qualitatively-observed wavelength at the time of this 
run; this was, as previously noted, an incident wave 


frequency cf 2.2 rad/sec. 
Тһе linear model replicates the cbserved data in 


pitch almost precisely, with the exception of the two nodes 


at 2.2 and 6.28 rad/sec, due to the dominant local wave 
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amplitudes in the observed data. 


As may be expected, pitch rate data is somewhat 
more noisy than pitch data, although not excessively so. 
ckserved pitch rate shows an unusual double peak in 
resonance, which also occurs їп the other two 
sea-conditions. The first peak occurs at 4.35 rad/sec, 
slightly higher in frequency than the observed oitch data, 
while the second peak occurs at 5.25 rad/sec and is of 
unexplained origin. This peak may be indicative cf the 
onset of the non-iinear behavior discussed in the preceding 


chapter, but that is only a speculative observation. 


The linear system pitch rate data shcws a 
Sharper rescnance than did linear pitch, the former being 
markedly peaked at 4.26 rad/sec, while the latter was 
relatively flat from 3.9 to 4.26 rad/sec. The agreement in 
the Sharpness of pitch rate response between the observed 
data and the linear system is quite good. Resonant rise 
also occurs earlier in the linear model, at approximately 
2.0 rad/sec; but there is no indication of the peak at 5.25 


rad/sec which occurs in the observed data. 


No significant nodes or nulls occur in this 


frequency range; the pitch node at n=1( w=3.15 rad/sec) does 
i 


not occur until an encounter frequency of 9.75 rad/sec. 


Overall, the observed and linear data ссарагг 


remarkably well under these conditions. 
b. Pitch and Pitch Rate-Cross Seas Run 


Despite the fact that  cross-seas runs induce 


modes of motion in roll and yaw that are unaccounted for in 
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this linear model, it was considered of interest to сспраге 
observed cross-seas data, Figs 40 and 41, with Figs 42 and 


43, the linear data for w=w. 
e 1 


The observed pitch FFT, Fig 40, is more 
distinctly resonant at 4.0 rad/sec under these sea 
conditions than in the ahead-sea run;response decreases 
continuously from the very low frequency end until almost 
the onset cf resonance at approximately 3.5 rad/sec. 
Resonance occurs at 4.0 rad/sec, and a distinct null cccurs 
at 5.25 rad/sec. The observed predominant wavelength, which 
as noted previously corresponds to an incident wave 
Irequency cf 2.2 rad/sec, is distinctly visible at that 


frequency. 


Coaparison with the linear system pitch response 
curves of Fig 42 is as striking as the ahead-seas runs, tne 
more so because of the modes of motion induced in the craft 
which are not in the linear model. AS in the FFT data, 
linear pitch response decreases until almost the onset of 
resonance, with a pronounced null at 3.5 rad/sec. This null 
does not appear in the observed pitch data, but is 
distinctly visible in the observed pitch rate data. Linear 
System resonance is sharper under these sea conditions, as 
is the cpserved data, and occurs at 4.26 rad/sec. Two 
high-frequency nulls at 5.25 and 7.0 rad/sec, and a node at 
5.0 rad/sec , are distinctly visible in the observed data, 


although less prominently. 


The observed pitch rate data shows a sharper 
rcll-off from the low frequency end of the spectrum, апа а 
flatter response thereafter untii the onset of resonance, 
than does pitch. The distinct doubie peak previously noted 
in the ahead-sea run, 15 distinctly visible in this run 


also. 
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Cemparısen ОЕ the pitch rate FFT, Fig 41, with 
the linear system pitch rate response of Fig 43 shows the 
same good agreement as did the pitch data; the low frequency 
roll-off characteristics are very much the same, and the 
null predicted at 3.5 rad/sec by the linear model is 
prcminently  visibie in the FFT. However, the higher 


frequency nodes and nulis are much less distinct. 
с. Pitch and Pitch Rate-Following Seas 


Due to the difficulties previously experienced 
in interpreting results of following seas runs with Bcoth's 
work, only a most general comparison of observed pitch and 
pitch rate FFT's, shown in rigs 44 and 45, and the linear 
system respcnses, shcwn in Figs 46 and 47, will be uade. 
Note that due to the technigue utilized in generating 
encounter frequencies in the linear model, the range of 
frequencies plotted in the linear model is much greater, (.3 
rad/sec to -40 rad/sec), resulting іп а much ccarser 
frequency resclution in the range of interest. Negative 
encounter frequencies, as noted in the preceding chapter, 


represent a modified ahead-sea condition. 


The linear system predicts a sharp roll-crf in 
pitch response from the very low frequency end to about -.25 
rad/sec, and a relatively flat response from there until tae 
onset of resonance, at about -2.0 rad/sec. Resonance cccurs 
quite sharply at -4.11 rad/sec. No nodes or nulls in this 
frequency range are predicted by the linear modei. The 
observed pitch data is in general agreement with this; in 
addition the predominant wavelength, which occurs under 
these sea conditions at an encounter frequency of -1.8 


rad/sec, may Бе observed in the FFT. 


The linear system and observed data do not agree 
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well in pitch rate, the linear system predicting an 
increasing response from the low-frequency end until 
resonance, while cbserved data shows, as should be expected, 
a low-frequency roll-off to a relatively flat mid-range 
response. No plausible explanation may be advanced for this 
variance at this time; however, considerable difficulty has 
been observed by beth Booth and this author in correctly 


interpreting following seas results. 


The observed pitch rate FFT shows, as have the 
runs in the cther two sea conditions, the mysterious double 


peak in pitch rate resonance. 


D. SUMMARY CF FREQUENCY DOMAIN ANALISIS 


Compariscn cf measured and linear system theoretical 
frequency response shows significant correlation under 
small-sea conditions; nulls predicted by the linear system 
may be readily observed in test craft data, and the 
frequency dcmain is relatively uncluttered by noisé. The 
bulk of noise in the system is not "white" noise, but tends 
to cluster about sharply-defined frequencies. The 
techniques developed in this chapter show the value of the 
frequency dcnrain as an interesting alternative to the time 


domain as a means of model verification and validation. 
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Figure 35 - REPRESENTATIVE WIDE-BAND NOISE SPECTRUM OF THE 
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Figure 36 - OBSERVED PITCH FREQUENCY SPECTRUM, 
15 KTS, AHEAD SEAS 
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Figure 37 - OBSERVED PITCH RATE FREQUENCY SPECTRUM, 
IS ETS AREAD SEAS 
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Figure 47 - LINEAR SYSTEM PITCH RATE FREQUENCY SPECTRUM, 
15 KIS, FOLLOWING SEAS 
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VII. SUMMARY 


CONE LUSTONS AND OBSERVATIONS 


The linear two-degree-of-freedom model developed in this 
work has been shown to be a reasonable approximation to both 
the Loads and Моїіопѕ ргоагап and to actual observed craft 
behavior, in both the time and frequency domain; its 
relatively uncluttered Simplicity has disclosed some 
interesting aspects of testcraft behavior, including the 
nulls and nodes that characterize the craft frequency 
response. Furthermore, this Simulation is extremely 
economical of CPU time, requiring less than three minutes of 
CPU time to generate the frequency response curves fcr up to 


100 discrete frequencies. 


The non-linear behavior described in Chapter V is 
somewhat supported by observations of testcraft data in the 
frequency domain in Chapter VI. This interesting behavior 


Mm apparently akin tio that ofa relaxation oscillator. 


Noise, which has interfered with varying degrees of 
severity in almost all transient analysis done with the 
recorded testcraft data at NPS, is reduced to little 
Significance using the frequency domain method applied for 
testcraft data reduction in this work. The major source of 
noise in tne XR-3 has been identified as primarily ignition 
noise from the eight gasoline engines aboard the craft, 


transmitted electrically to the data acquisition systen. 
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В. RECOMMENLATIONS FOR FUTURE STUDY 


The linear heave-pitch model developed in this work can 
readily be expanded to include roll, yaw and speed effects. 
Completicn of this expansion would result in an analytical 
Six-degree-cf-freedom model, and lead to a more thcrough 
understanding of CAB craft behavior. Ёп example of the 
potential of this approach may be seen in ref [9], in which 
a rudimentary pitch control system, based upon this model 
and utilizing optimal Ricatti solutions, promises to reduce 
sea-state pitch accelerations to less than ten percent of 
the open lcop response. The practical problems involved in 
the isplementation of a test-bed ride control system fcr the 
XR-3 testcraft do not appear to be severe, aside from tne 
physical installation of a suitable set of stern-plane 
control surfaces; an observer, based ona linear model, 
could be utilized as an optimal observer in a Kalman 
filtering scheme to generate some of the system states which 
may not otherwise be suitable for feedback. Since TEAS 
system is only piecewise linear, the microprocessor appears 
as the most feasible means of implementing this control 
system, storing the optimal gains for several operating 
speeds in a table look-up or interpolation method. Most of 
the systen development could be implemented оп а 
slightly-mcdified version of the Loads and Motions prcgran. 
This avenue has great future potential for expanding the 
sea-state envelope and improving the habitability of САВ 


Craft. 


The linear freguency domain response simulation should 
prove useful in the study of the non-linear  betavior 
described in Chapter V. It is felt that this non-linearity 
can be modelled in the complex plane utilizing a describing 
function.  Fcrtunately, the complex gains of the linear 


system are directly available from the linear system in the 
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program developed. This is an extremely interesting area of 


future study. 


In the area of parametric analysis, such a linear model 
may be utilized with the technique of root locus analysis. 
An elementary program was developed for this purpose and is 
included at the end of the text with other programs 
developed for this investigation. However, time and space 
did not permit its inclusion into the main text. 
Furthermore, such parametric analysis was inconsistent with 
the main thrust of this effort. It is nevertheless hoped 


that the program will be of benefit in future studies. 


The XR-3 testcraft was not available for data collection 
throughout the spring of 1977. The frequency domain data 
reduction techniques developed in this work could not, 
therefore, be utilized in anything approaching an in-depth 
study of craft behavior. One area in which these frequency 
domain techniques may prove most useful would be the 
investigation of actual testcraft gains in heave at various 
frequencies, to ascertain the significance of wave-plenun 
Chamber interactions. The 15 dB discrepancy between the 
Loads and Motions program, which assumes a plenum wave, and 
this linear model, which does not, indicates this to bea 


vital factor in ths accuracy of any model. 


A serious effort at ignition filtering is certainly in 
order. Ignition suppression and power-line filtering is a 
relatively straightforward approach that can greatly enhance 
the quality of data collected from the XR-3. 


An additional technique for data reduction to the 
frequency domain may become available in the ver near 
Future. The MDS-8 microprocessor-based data acguisition 
system, developed by Englehart in ref рој, will be capable 


of sampling sequentially and storing in mass memory up to 14 
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analog channels, with a sampling rate up to 7000 samples per 
second. Upon completion of the sampling and storage, the 
MDS-8 can perform a variety of data manipulations upon the 
data file, including an FFT; the data file, both in the raw 
form and after manipulation, is available to the NPS IBM 
360/67 computer installation via the CP/CMS time-sharing 
link. The possibilities of this system as a data reduction 
system are, compared with technigues previously utilized DOM 
data reducticn, impressive. However, much remains to be 


developed in this area. 
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APPENDIX A 


NIDEBAND NOISE DATA 


This appendix is a resume of the more significant 
wide-band ncise spectra observed in the FFT's of XR-3 
testcraft data collected at Lake San Antonio, 20 January 
1977. 


The following plots are FFT's of testcraft data, reduced 
as described in Chapter V, for a maximum upper frequency of 
200 Hz, and a frequency resolution of 512 samples. 
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This program provides a root locus solution tc the 
two-degree-of-freedom linear model. The reserved СМР 
variable TIME is redesignated as the parameter of interest 
with the RENAME statement, with the range of variation 
specified on the FINTIM card. One additional variable may be 
stepped, using the PARAMETER statement. 


The parametric eigenvalue program also must access the 
IMSL Subroutine Library; see the Frequency Domain computer 
pregram listing for the appropriate ӘСІ tc acheive this 
linkage at tke W.R. Church Computer Facility. 
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